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Abstract. We present ground-based 2.1 to 4.2 /im observations of four low-mass protostars. We searched for the 
4.1 ^m OD stretch band, characteristic of solid HDO in grain mantles. We did not detect solid HDO in any 
of the four sources, but we derive 3a upper limits from 0.5% to 2% for the HDO/H2O ratio depending on the 
source. These ratios provide strong constraints to solid-state deuteration models when compared to deuterium 
fractionation values observed in the gas phase. We discuss various scenarios that could lead to such a low water 
deuteration compared to the high formaldehyde and methanol deuteration observed in the gas-phase. 
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1. Introduction 

Large deuteration fractionations — ratio of deuterated 
over normal isotope — are a common characteristic of 
dark, dense molecular clouds, with abundances of deuter- 
ated species some 0.01-0.1 of the normal isotopomer (e.g. 
ITurnerl 2001). This deuteration is thought to reflect the 
small zero point vibrational energy differences between a 
deuterated species and its fully hydrogenated counterpart, 
which drives the chemistry towards deuterated species at 
low temperatures (10 K: IWatsonI 1973). In view of their 
elevated temperatures (> 100 K), the warm gas around 
protostars shows unexpected high abundances of deuter- 
ated species. Fractionations of lO^'^ to more than 0.1 
have previously been observed. Low-mass protostars seem 
to show even greater fractionations than high-mass pro- 
tostars. Indeed, observations of the low-mass protostar 
IRAS 16293-242 2 showed that D2CO / H2CO - 10 % 
IjCeccarelli et al.l 1998), a fractionation about 25 times 
larger than in the Hot Core in Orion l|Tiirnerl 1990). 
Similarly large amounts of doubly-deuterated formalde- 
hyde have subsequently been observed towards a sam- 
ple of low-mass protostars IjLoinard et al.l 2002). Singly- 
and doubly-deuterated methanol were detected towards 
IRAS16293-2422 UParise et al.l 2002), showing fraction- 
ation ratios of 0.9 (CH2DOH) and 0.2 (CHD2OH). 
Deuterated methanol was also observed towards a sample 
of other low-mass protostars IjParise et al.l in prep) , con- 
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firming the huge deuteration in these objects, whereas the 
fractionation ratio CH2DOH/CH3OII is only 0.04 in the 
Orion high-mass star-forming region Pacq et aT] 1993). 

Because the high temperatures of the gas around pro- 
tostars prevent significant deuterium fractionation, the 
observed high deuterations are thought to reflect a pre- 
vious cold phase. Molecules that formed during the dark 
cloud phase — either in the gas phase or on the grain 
surfaces — are believed to be stored in an ice mantle 
on dust grains, which evaporates once the YSO heats 
its environment above the ice sublimation temperature 
l|Ceccarelli et al. 1 2001). Presumably the current gas phase 
has not had time to return to equilibrium, which takes 
some 3.10"^ years | |Charnley et al.| 1992, iCaselli et al. 1 1993. 
ICharnley et"aLl l997). 

A previous detection of solid-state HDO in the high- 
mass luminous W33A and NGC7538 IRS9 protostars has 
been claimed by Teixe ira et al.l (1999). However, recent 
analysis of the same data seems to challenge this detec- 
tion, giving an upper limit of 10~^ for solid HDO/H2O 
in NGC7538 IRS9 (Dartois '^ran 2003V Subsequent VLT 
observations of W33A also give an upper limit for the 
deuterium fractionation of ~ 10~^ (jDartois et al.l 2003). 
a value still consistent with the typical deuteration ob- 
served for the gas phase of hot cores around high mass 
protostars. The search for HDO in high-mass protostars 
was thus thought to be inconclusive largely because the 
expected fractionation was very small. Because low-mass 
YSOs show gas-phase deuterium fractionations almost 100 
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times higher in their hot cores, the sohd HDO feature is 
expected to be much stronger, and motivated this study. 
Presently published grain chemistry models predict an 
HDO fractionation of about 40% (i.e. almost all the deu- 
terium would be locked in HDO) if deuterated water forms 
on the grain surfaces at the same time as formaldehyde 
and methanol l|Caselh et al.l 2002, IStantcheva & HerbstI 
2003). 

In this paper, we present a search for the HDO stretch- 
band at 4.1 /im towards four low-mass protostars. In sec- 
tion 2 we describe the observations, and derive the H2O 
ice column density and upper limits for HDO column den- 
sities in section 3. We discuss the derived fractionation 
ratios, comparing to available gas-phase observations (for 
SVS12 and SVS13) in section 4, and conclude in section 
5. 

2. Observations 

2.1. Source selection 

As noted in the Introduction, the previous searches for 
solid HDO have been carried out towards high-mass pro- 
tostars, with very modest success. In this study we focused 
on low-mass protostars, where gas-phase deuteration has 
been observed to be much higher than in high-mass pro- 
tostars. 

The sources of the present study were selected to be 
bright enough at NIR wavelengths, so that high enough 
S/N can be obtained to detect weak absorption features 
against the continuum, and to have large J-K colour index, 
which would indicate a high extinction. 

The first criterion, by definition, excludes the most em- 
bedded low-mass protostars, the Class sources, because 
the NIR continuum is too faint, leaving the more evolved 
Class I sources. Even though not totally conclusive, the 
study of ILoinard et al. l (2002) on a sample of four Class 
I sources seems to indicate that also in those sources the 
molecular deuteration, and specifically the D2CO/H2CO 
ratio, remains relatively large and comparable to that 
found in Class sources. Encouraged by this, we selected 
the four sources in Tabled Two of our sources were ob- 
served bv lLoinard et al.l f2002). In the following we give a 
brief description of each selected source. 

NGC1333 SVS12, also called IRAS6, is proba- 
bly a Class I source of ^ 28 Lq in the IRAS bands 
Pennings et al.| 1987) assuming a distance of 350 pc. New 
distance estimates of the NGC1333 complex tend to put 
it closer at about 220 pc l|Cernisl l990). which would bring 
the SVS12 luminosity to ~ 10 Lq. The infrared source is 
close to the Herbig Haro object HH12, although SVS12 
probably is not the exciting source of the HH12 flow. 
Millimetre and submillimetre maps resolved the source 
in possibly three components IjSandell fc Kneel 2001). A 
detailed CO 2-1 map of the region did not detect any out- 
flow emanating from this source fKnee fc S andcU 2000) so 
it may be a relatively evolved source indeed, and possibly 
a background Class II source. In this case, we may be ob- 



serving the deuteration in the cloud ices rather than in 
the circumstellar material. ILoinard et al.l (2002) detected 
abundant doubly deuterated formaldehyde, D2CO/H2CO 

NGC1333 SVS13, also called IRAS03259-I-3105 
dJennings et al.| 1987), is the best studied source of the 
region, thanks to its well coUimated flow of HH objects, 
the HH7-11 complex. It is a variable object and one of the 
brightest NIR sources of the region, having an average 
bolometric luminosity of 115 Lq IjMolinari "6^0 1993: 45 
Lq if the distance is 220 pc), and it is very probably a mul- 
tiple system (|Lefloch et al.l 1998a. b: .Bachiller et all 1998: 
Loon ey et al.| 2000). Also in this source ILoinard et al.l 
(2002) detected abundant doubly deuterated formalde- 
hyde, D2CO/H2CO ~4%. 

L1489 IRS (IRAS04016-h2610) is a low-luminosity 
source (Lf,oi=3.7 Lq), which acquired recent attention 
because it is suspected to be one of the rare cases 
of a source in transition from Class I to Class II 
CHo gerheijde fc Sandell| 2000; 'Hoge rheijde] 2001). It ap- 
pears to be surrounded by a relatively massive and young 
disk ( |Boogert et al.| 2002a), whereas the envelope seems 
to be largely swept out. 

TMRl (IRAS0436H-2547) is a prototypical Class I, 
low-luminosity source (Lbo;=2.9 Lq: Hogcrhc ijde et al.| 
1998). Both L1489 IRS and TMRl belmig to the Taurus 
complex at a distance of 140 pc. 

2.2. Observations 

Using the SpeX instrument ( |Rayner et al.| 2003) on the 
IRTF telescope on Mauna Kea (Hawaii), we observed the 
four low-mass class I protostars listed in Table The ob- 
servations were performed in December 2002 under condi- 
tions of good seeing (0.6" at K, on average). We used the 
2.1-5.0 iim cross-dispersed mode, which acquires the full 
spectrum between 2.1 and 5.0 /xm simultaneously. We used 
the 0.3" slit, which provides a spectral resolution of 2500, 
for the first three objects, and the 0.5" slit (R=:1500) for 
TMRl. Standard stars were chosen from the IRTF AOV- 
star database to be as bright and as close to the studied 
object as possible, and are listed in Tabled The airmass 
difference between the object and the associated standard 
was always less than 0.1. 

The observations were reduced using the Spextool 
software IjVacca et al.l 2003, Cushing et al. in prep), 
which is available on the IRTF SpeX webpage 
I http://irtfweb.ifa.hawaii.edu/Facility/spex/ll. The full 
spectra for the four protostars, smoothed to a resolution 
of ~ 10 cm~^, are presented in Fig.Q] 

The HDO absorption feature was not detected in 
any of the four sources. We derive upper limits for 
the HDO/H2O ratio in next section. The spectra reveal 
two hydrogen lines at 2468 cm^^ (4.05/im) and 2673 
cm~^ (3.74^m), as well as CO emission lines (SVS13, 
ICarr et aP 1992). The deep absorption feature at 3300 
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HD20995 


NGC1333 / SVS13 
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15 min 
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9.3 
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20 min 


HD23258 
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0.08 
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10.54 


0.5" / 1500 


53 min 
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Table 1. Observed sources, coordinates, NIR magnitude, J-K colour index (indicative of extinction, cf text), K 
magnitude and observation parameters. 
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Fig. 1. Spectra of the four protostars, smoothed to 10 cm ^ resolution. The vertical lines delineate the region of strong 
atmospheric absorption features. 



cm ^ (3.03/im) is the OH stretch band. On its right wing, 
the so-called "3.47 /xm feature" (2882 cm^^) is visible. 

3. Water features 

3.1. The OH stretch band 

All four spectra show a prominent OH absorption feature 
at 3300 cm^^ (3.03 fim). The abundance of sohd hydro- 
genated water can be derived from this OH stretch band. 
We first define the continuum Icont around the OH feature. 
We can then plot the optical depth r defined by : 



The column density of H2O molecules can be derived 
from the optical depth with the following relation : 



Nh,o = 



A 



H20 



where the band strength AH2O = 2 x 10 cm/molecule 
l|Dartois| 1998) and / rdi' is expressed in cm^^. 

3.1.1. SVS12, L1489 IRS and TMRl 

For these three sources, we defined the local continuum 
around the OH feature by fitting a second order polyno- 
mial to the frequency ranges [2400-2500, 2700] and [4000, 
4400-4600] cm~^. These frequency ranges were chosen to 
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Fig. 2. Example of determination of the continuum for 
the study of the 3.0/xm feature. The dotted lines show the 
frequency ranges used to fit the continuum 

avoid the region around the 2882 cm~^ (3.47 iim) feature 
in the continuum fit. An example of the derived continuum 
is shown in Fig. 13 

The water feature for the three sources shows no struc- 
ture and is thus attributable to amorphous ice. We in- 
tegrated the contribution of the band between 3000 and 
3555 cm^^. This integrated optical depth corresponds to 
a lower limit on the H2O column density. The OH bands 
for these three sources are plot in an optical depth scale 
on Fig. 13 

The optical depth for the 3300 cm^^ feature for each 
source and the derived H2O column density are summa- 
rized in Table H Our optical depth for L1489 IRS IS con- 
sistent with the observation at low resolution (R=150) by 
ISato et al.i ri990: t = 2.9 ± 0.2), but we find a lowe r (2.0) 
optical depth for SVS12 (r = 2.9 ± 0.2. ISato et aTl 1990V 
This inconsistency certainly arises from the continuum de- 
termination, as lSato et al.l onlv had observations between 
2.4 ^m (4167 cm"!) and 3.8 ^m (2632 cm'^). 

3.1.2. SVS13 

SVS13 appears to be a more unusual source than the 
other three. First of all, its OH stretch band has a non- 
gaussian shape, which suggests that the ice is in a crys- 
talline state. Second, the spectrum is characterized by a 
strange inflexion in the frequency range [3900, 4700] as 
well as [2400, 2600] cm'^ The inflexion in the [2400, 2600] 
cm^^ range cannot come from the standard star, as we 
used the same standard object for SVS12, which shows no 
inflexion. Moreover, ISO-SWS observations at low resolu- 
tion of this source (ISO archive tdt 65201959) have the 
same inflexion. This strange behaviour of the spectrum 
makes it difficult to define the local continuum. 

We defined several different continua, by fitting first or- 
der polynomials on different anchoring frequency ranges. 
An example is shown in Fig. 2J In all cases, the optical 
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Fig. 3. OH stretch band for SVS12, L1489 IRS and 
TMRl. For clarity, the plots for TMRl and L1489 IRS 
have been offset by 0.4 and 0.7 respectively. The verti- 
cal lines delineate the band where the integrated optical 
depth was computed. 



depth plots have a broad absorption on the right side of 
the OH feature. This absorption could be due to scatter- 
ing in the OH wing IjDartois et al.1 2002'). The integrated 
optical depth varies from 190 to 220 cm^^. We chose the 
smallest value for a lower limit on the H2O column density. 

3.2. Upper limits on the HDO abundance 

To determine upper limits on the HDO abundances, 
we first subtracted the two hydrogen lines at 2468 and 
2673 cm^^. To do so, we fitted a first or second order 
polynomial in the [2400, 2800] cm~^ frequency range, ig- 
noring a small frequency range around each of the two 
hydrogen lines (horizontal marks in Figure[5l). Continuum 
fits are plotted for each source in Figure [3 The hydro- 
gen lines were subtracted from the data by replacing the 
data in the associated frequency ranges by the interpo- 
lated continuum. 

We then defined a local continuum around the ex- 
pected 2457 cm~^ (4.07/xm) OD stretch band by fitting 
a first or second order polynomial in the [2410,2440] and 
[2500,2800] cm~^ frequency range. 

An upper limit on the HDO column density was de- 
rived taking into account that the expected feature is 
rather broad (120 cm~^). The data were smoothed ac- 
cordingly to a resolution of 30 cm^^. The upper limit on 
the HDO column density was then derived by estimat- 
ing the biggest gaussian feature fitting in a 3a-wide band 
around the polynomial fit. 

The gaussian, centered on 2457 cm~^, was chosen to 
have a fixed FWHM of 120 cm~^ (0.2 /im), as measured in 
the laboratory for amorphous HDO ffiartois et al.l 2003). 
The extreme gaussians for each source are plotted on an 
optical depth scale in Figure [S] (see insets) . The dashed 
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Source 


TOH 


r Tn M dv 

(cm-') 


N(H20) 
(10^» cm-2) 


roD 


'(cm-i) 


N(HDO) 
(10''' cm-2) 


HDO/H2O 


NGC1333 / SVS12 


2.0 


> 750 


> 3.8 


< 0.005 


< 0.63 


< 1.8 


< 0.5 % 


NGC1333 / SVS13(a) 


0.55 


> 190 


> 1.0 


< 0.016 


< 2.0 


< 5.6 


< 5.9 % 


NGC1333 / SVS13(b) 








< 0.022 


< 0.58 


< 1.6 


< 1.7 % 


L1489 IRS 


3.0 


> 1160 


> 5.8 


< 0.013 


< 1.6 


< 4.4 


< 0.8 % 


TMRl 


2.0 


> 770 


> 3.9 


< 0.013 


< 1.6 


< 4.4 


< 1.1 % 



Table 2. OH and OD stretch band characteristics. Upper limits for HDO parameters are 3-cr. (a) in the amorphous 
case, (b) in the crystaUine case (see text). 
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Fig. 5. OD stretch band for the four protostars. The horizontal bars show the frequency ranges ignored in the con- 
tinuum fit. The continuum is drawn as a dashed line. In insets, the data are plotted on an optical depth scale versus 
frequency, after smoothing to 30 cm^^ resolution. The dashed lines are the 3ct limits on the flux, transposed on an 
optical depth scale. The gaussian shows the maximum amorphous HDO absorption. 



lines represent the ± 3(t limits associated with the contin- 
uum determination at 30 cm^^ resolution. 

The upper limit on the HDO column density is then 
derived from the relation : 

r Tciv 

Nhdo < 

AoD 

where / rdv is the area of the extreme gaussian, and Aqd 
was taken to be 3.6 x 10~^^ for a conservative upper 



limit estimate {Aqd — (4.3 ± 0.7) x 10 -^^cm/molecule, 
IDartois et all 2003). 

For SVS13, we estimated an upper limit on the HDO 
column density for the case of crystalline ice. The gaussian 
was then centered on 2427 cm~^ (4.12 /xm), with a FWHM 
of 30 cm~^ (0.05 fim) following laboratory experiments 
HDartois et al. 1 2003). In that case, the data were smoothed 
to 10 cm~^ resolution, the results for which are plotted on 
Figure E| 
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Fig. 4. Example of determination of the continuum for 
the study of the 3.0/xm feature. The dotted hues show the 
frequency ranges used to fit the continuum. 

The 3cr upper limits for the HDO column densities and 
associated HDO/H2O ratios are hsted in Table [3 

4. Discussion 

The observations presented here represent the first at- 
tempt to study the solid-phase deuteration in low-mass 
protostars, namely sources with luminosities of a few tens 
of solar luminosities. Previous studies focused on sources 
whose luminosity is larger than about 200 (e.g. Dartois 
et al. 2003), and where the deuteration of gas phase 
molecules is relatively low. Specifically, the gas-phase 
D2CO/H2CO ratio is < 0.4% in high-mass protostars, 
whereas it is > 4% in low-mass protostars (|Loinard et al.l 
2002). For this reason, the present observations represent 
the most stringent constraint and a crucial test on the re- 
lation between the solid H2O deuteration and the large 
deuteration observed in the minor gas-phase species. 

Loinard et al. (2002) reported the observation of 
D2CO/H2CO in two (SVS12 and SVS13) out of the 
four sources studied here. In both sources the measured 
D2CO/H2CO ratio is around 4%, a value similar to 
that found also in the other two Class I sources stud- 
ied by ILoinard et al.l (2002), and just a factor two or 
three smaller than the average value measured in Class 
sources. Hence the case of SVS 12 and SVS 13 is likely 
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Fig. 6. OD stretch band for SVS13. The horizontal bars 
show the frequency ranges ignored in the continuum fit. 
The vertical bar is a representative ± 3(t error bar. The 
maximum gaussian absorption, superimposed as a solid 
line, has been computed under the assumption of crys- 
talline ice (see text). 



representative of the formaldehyde deuteration in Class I 
sources, and more generally in embedded low-mass pro- 
tostars, within a factor two or three. Note that the 
HDCO/H2CO ratio in embedded protostars is even larger, 
~ 10%, as observed in the Class source IRAS16293- 
2422 IILoinard et al.l 2000. Ivan Dishoeck et al.1 1995'1. Not 
to mention the case of the deuterated methanol, which 
is about as abundant as the main isotope l|Parise et al.l 
2002). 

The present observations clearly show that in low-mass 
protostars the solid HDO/H2O ratio is lower than about 
2% (3-0-), definitively lower than the observed fractiona- 
tion of the gas-phase formaldehyde. This is a crucial step 
towards the understanding of the deuteration from an ob- 
servational point of view. 

Before discussing the consequences of the compari- 
son of these ice observations with gas-phase observations, 
it is important to understand where these ices are lo- 
cated along the line of sight and whether the ices we 
are probing have experienced the same chemical history 
as the ices that evaporated in the hot core and released 
the high-deuterated molecules such as formaldehyde and 
methanol into the gas-phase. In principle, a substantial 
part of the ice absorption may arise in the foreground qui- 
escent gas not directly associated with the protostar itself 
dBoogert et al.| 2002b). Observation of the pure sohd CO 
band at 4.67 /xm may help elucidate how much thermal 
processing the ices have undergone and hence their prox- 
imity to the YSO. Indeed, pure CO absorption will only 
be present for ices less than 15 K and hence is expected 
to arise only in the dark foreground cloud. The CO band 
towards L1489 IRS has been observed and modelled by 
[Boogert et al.| (2002a). They find an optical depth of 0.4 
for the pure CO component (see their figure 5). The distri- 



Parise et al.: Search for solid HDO in low-mass protostars 



7 



button of pure solid CO in the Taurus molecular cloud has 
been well studied (jWhittet et al.l l989. Chiar et al.l 19951. 
The optical depth in the unperturbed dark cloud scales 
linearly with the total extinction {Tpureco = 0.08(Ai, —6), 
KJhiar et al.l 1995). The pure CO optical depth towards 
L1489 IRS corresponds then to a visual extinction of 11 
for the foreground cloud and a visual extinction of 18 for 
the envelope associated with the protostar. The H2O ice 
optical depth in the Taurus dark cloud shows a similar lin- 
ear correlation with A„ {TH2O — 0.08{Ay — Sl. lChiar et al.l 
1995). Hence, the H2O optical depth associated with the 
foreground cloud is only 0.6 and ice grain in the envelop 
of L1489 IRS contribute TH20,ice — 2-4. Thus, our upper 
limit on the fractionation of H2O ice in the immediate 
environment of this protostar is 1%. 

In the following we discuss the consequences of 
these two sets of observations, gas-phase and solid-state 
deuteration. We start discussing two possibilities: either 
formaldehyde (and methanol) forms at the same time as 
water ices or water ices are formed in a previous phase and 
therefore they arc imprinted with the -different- deuter- 
ation degree of this phase. 

There is now a general consensus, supported by 
several observational and theoretical works, that both 
formaldehyde and methanol form on the grain sur- 
faces, probably after hydrogenation of CO molecules 
that have stuck onto the grains (e.g. theoretical stud- 
ies of IShalabie a fc Greenb'eri 1994; ITielens & Whittetl 
19 97; [Charnley et al.| 1 997, and laboratory experiment 
of IWatana^^^Couchil 2002). although one laboratory 
experiment seems to challenge this last hypothesis 
(Hiraok a et al.l 2002^1. If we assume now that water ices 
are formed at the same time by hydrogenation of oxy- 
gen atoms sticking on the grain surfaces, the deutera- 
tion degree is set by the atomic D/H ratio of the accret- 
ing gas (e.g. ITielensI 1983; ICaselli et al.l 2002), and can- 
not be substantially different in the three species (water, 
formaldehyde and methanol), unless selective deuteration 
is at work in formaldehyde and methanol with respect to 
water. Detailed chemical models have been developed for 
the deuterium fractionation on grain surfaces. These lead 
to enhanced fractionation of formaldehyde with respect to 
that in other ice species centering on H-abstraction reac- 
tions, driven by the zero-point energy difference (Tielens 
1983). However, our insights in grain surface chemistry 
have evolved considerably since then and these models 
have not been fully revisited. 

In the second hypothesis, water, formaldehyde and 
methanol would be formed during different phases char- 
acterized by different physical conditions which then may 
lead to different levels of fractionation (cf.. IDartois et aO 
2003). In particular, if the methanol and formaldehyde 
are formed during the later stages of accretion - when 
most of the gas has accreted already -, their fractiona- 
tion could be high IjR.oberts et al.l 2003). Such a model 
would require that the physical conditions in the accret- 
ing gas which favor deuteration also favor formaldehyde 
and methanol formation but not water formation. This 



hypothesis encounters some problems as the high abun- 
dance of formaldehyde and methanol in hot cores is gen- 
erally thought to reflect the evaporation of water-rich ices 
HLoinard e t al.' 2000; ' Ceccarelh etaTl 2001, ISchoier et all 
2002). Nevertheless, observation of abundant D2CO in 
the outer envelop of IRAS16293 l|Ceccarem et al.l 2001) 
provides some observational support for onion-like man- 
tle structure, with deuterated species trapped in the CO- 
rich ices that evaporate around 15K. In this scheme gas- 
phase observation of deuterated molecules would trace the 
deuteration ratio in these CO-rich ices, whereas the solid- 
phase observation would trace the bulk of the H20-ices, 
where deuteration may be less important. 

Alternatively, water ices are formed by condensation 
of water, copiously formed in molecular shocks, occuring 
during the cloud phase ( |Bergin et al.| 1999). Also in this 
case the HDO/H2O ratio would be low, reflecting the high 
temperatures (~ 300 K) at which water would have been 
formed. Within this scenario, formaldehyde and methanol 
are not shock-produced and hence show quite different 
(higher) deuterium fractionations characteristic of grain 
surface (or gas phase) formation at low temperatures dur- 
ing the preshock phase. Of course, shocks also have other 
chemical consequences - notably enhanced abundances of 
SiO and SO are considered to be common tracers of shocks 
(e.g., Codclla et al. 2002; Schilke et al. 1997, and refer- 
ences therein). However, this pertains to gas phase species 
and, moreover, cannot address the deuteration issue at 
hand here. 

None of these models seem to be fully developped and 
all seem to have some difficulty explaining the observa- 
tions. Perhaps, it is time to reconsider models for the se- 
lective deuteration of formaldehyde and methanol on grain 
surfaces. 

5. Conclusions 

We presented a search for solid HDO in grain mantles to- 
wards low-mass protostars. We did not detect HDO but 
derive upper limits on the HDO/H2O frationation ratio 
of 0.5 to 2%. These upper limits definitely show that 
solid water is much less deuterated than other molecules 
observed in the gas phase, such as formaldehyde and 
methanol. The origin of these differences is not fully un- 
derstood. 
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